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ABSTRACT 

Gas within molecular clouds (MCs) is turbulent and unevenly distributed. Interstellar shocks 
such as those driven by strong fluxes of ionising radiation (IR) profoundly affect MCs. While 
small dense MCs exposed to a strong flux of IR have been shown to implode due to radiation- 
driven shocks, a phenomenon called radiation driven implosion, larger MCs, however, are 
likely to survive this flux which in fact, may produce new star-forming sites within these 
clouds. Here we examine this hypothesis using the Smoothed Particle Hydrodynamics (SPH) 
algorithm coupled with a ray-tracing scheme that calculates the position of the ionisation-front 
at each timestep. We present results from simulations performed for three choices of iR-flux 
spanning the range of fluxes emitted by a typical B-type star to a cluster of OB-type stars. The 
extent of photo-ablation, of course, depends on the strength of the incident flux and a strong 
flux of IR severely ablates a MC. Consequently, the first star-formation sites appear in the dense 
shocked layer along the edges of the irradiated cloud. Radiation-induced turbulence readily 
generates dense filamentary structure within the photo-ablated cloud although several new 
star-forming sites also appear in some of the densest regions at the junctions of these filaments. 
Prevalent physical conditions within a MC play a crucial role in determining the mode, i.e., 
filamentary as compared to isolated pockets, of star-formation, the timescale on which stars 
form and the distribution of stellar masses. The probability density functions (PDFs) derived 
for irradiated clouds in this study are intriguing due to their resemblance with those presented 
in a recent census of irradiated MCs. Furthermore, irrespective of the nature of turbulence, the 
protostellar mass-functions(MFs) derived in this study follow a power-law distribution. When 
turbulence within the cloud is driven by a relatively strong flux of IR such as that emitted by 
a massive O-type star or a cluster of such stars, the MF approaches the canonical form due to 
Salpeter, and even turns-over for protostellar masses smaller than ~Q.2 M . 

Key words: Ionising radiation - Hllregions - Filaments - Star-formation - Initial mass func- 
tion 



1 INTRODUCTION 

Stars condense out of cold dense clumps of molecular gas called 
prestellar cores which appear to be located in much larger clouds, 
often filamentary in shape. While the formation of stars themselves 
is relatively well understood, a number of questions such as the 
distribution of stellar masses or even the formation of prestellar 
cores remain amongst the unsolved questions in contemporary as- 
tronomy. It was commonly believed that stars belonged to larger 
clusters populated with several tens of other members (e.g. Palla 
2005), but the belief is now contested by more recent observations 
of dozens of clusters (e.g. Guttermuth et al. 2009). 

A better understanding of the processes leading to the as- 
sembly of potential star-forming gas probably holds the key 
towards unravelling issues fundamental to the theory of star- 
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formation. Primarily, they are related to the mode of star-formation, 
in other words clustered as against isolated or a mixture of the two. 
The origin of the stellar IMF, its universal nature and the possible 
relation with the distribution of core masses is another issue that 
has baffled researchers. Also a related question is that about the 
time-scale on which a MC is likely to form stars. For a more de- 
tailed discussion of these points the interested reader is referred to 
Hartmann et al. (201 1) and Bonnell et al. (201 1). 

Furthermore, observations in a number of wave-bands, and in 
particular, the molecular lines at sub mm-wavelengths have now 
demonstrated the turbulent interiors and the non-uniform distribu- 
tion of gas within molecular clouds (MCs). In fact, a significant 
fraction of molecular gas appears to reside in dense filamentary 
regions (e.g., Elmegreen 1997, Nutter et al. 2008, Andre et al. 
2010). However, numerical simulations over the last decade have 
demonstrated with considerable success the tendency of turbulent 
gas to generate filamentary structure on a relatively short times- 



© 2002 RAS 



2 Anathpindika. S. and Bhatt. H. C. 



cale, usually smaller than a sound-crossing time (e.g., Klessen et 
al. 2000, Jappsen et al. 2005 amongst a number of other works). 
Though turbulence provides support against self-gravity on a global 
scale, locally, on scales comparable to the length of driven modes, 
it assembles gas that can eventually become self-gravitating which 
probably explains star-formation within larger filamentary struc- 
tures. Young protostars compete with each-other to acquire mass 
from their natal pool, a scenario better-known as competitive accre- 
tion (Bonnell et al, 2001). Evidently the formation of these dense 
filaments and their possible fragmentation depends on the prevalent 
physical conditions within the parent MC. 

Numerical simulations by Ballesteros-Paredes et al, (1999), 
Heitsch et al, (2008), and Anathpindika (2009 a,b) amongst those 
by several other authors, have shown that dense filaments of gas 
form rapidly via fragmentation of larger pressure-confined gas bod- 
ies that may form due to collisions between turbulent flows. The 
fragmentation itself is a result of a process called gravoturbulent 
fragmentation, an interplay between the gravitational instability 
and hydrodynamic instabilities. These simulations have shown that 
the star-formation process, beginning with the assembly of gas in 
dense pockets, once triggered, proceeds rapidly due to the non- 
linear growth of instabilities. However, the fact that a MC prob- 
ably spends considerable time in a relatively quiescent state before 
star-formation is triggered, according to some authors supports the 
hypothesis of slow star-formation (e.g., Krumholz & Tan 2007). 

MCs in astrophysical environments though, are also exposed 
to ionising radiation (IR), which is likely to profoundly affect 
their evolution and possibly trigger new episodes of star-formation. 
Large nebulae such as the Rossette nebula (e.g. White et al. 1997; 
Poulton et al. 2008), the Eagle nebula with its famous pillars of 
creation and trunks (e.g. Sugitani et al. 2002), the M17 nebula (e.g. 
Hoffmeister etal. 2008), the star-forming cloud NGC7538 (e.g. Ojha 
et al. 2004), the nebula SF038 in IC1396 (Choudhary et al. 2010), 
and cometary globules such as IC1848 (e.g. Lefloch et al. 1997), and 
CG12 (Haikala & Reipurth 2010), are just a few examples of star- 
formation triggered by the IR. The Rosette MC, for instance, is irra- 
diated by the nearby star-cluster NGC2244 and the radiation induced 
turbulence appears to have generated a network of filaments within 
the main cloud. More recent observations of the Rosette MC using 
the multi-waveband infrared cameras, PACS and SPIRE, on-board 
the Herschel space-observatory have revealed young star-forming 
regions in the junctions of these filaments (Schneider et al. 2012). 

Secondary star-formation, i.e., cases of new episodes of star- 
formation triggered by the flux of IR emitted by an earlier gen- 
eration of stars have also been reported in for e.g., the Orion MC 
(Wilson et al. 2005). Another prognosis, observed sometimes in 
regions of triggered star-formation, is the increasing age of stel- 
lar population along a certain direction, also called the sequential 
mode of star-formation (e.g., Maaskant et al. 2011). Hot radiation 
from young stars is therefore believed to play a crucial role in con- 
trolling the rate of galactic star-formation, and furthermore, it is 
also an important source of turbulence in galactic disks (e.g. An- 
drews & Thompson 201 1 ; Krumholz & Matzner 2009). 

Over the last three decades several authors have studied ana- 
lytically the likely fate of MCs irradiated by a strong flux of IR. 
Dyson (1973), for example, suggested flattening of an irradiated 
cloud, whereas according to Bertoldi (1989) and Bertoldi & McKee 
(1990), an irradiated cloud was more likely to implode under the 
influence of a radiation driven shock-wave within the cloud, a phe- 
nomenon known as radiation driven implosion. The cloud so ac- 
celerated was shown to acquire a cometary structure with a charac- 
teristic velocity of the order of a few km/s. Stars could form within 



this comet-shaped globule (CG) . Indeed, a number of such globules 
have been observationally reported as cited above. The formation 
of CGs due to the RDI has been demonstrated numerically for small 
dense clouds by for instance, Lefloch & Lazareff (1994), and more 
recently by Bisbas et al. (201 1). 

Dale et al. (2007), on the other hand, showed that the expos- 
ure of a massive self-gravitating, turbulent MC to a flux of strong 
IR significantly reduced the gas-to-star conversion efficiency. While 
these authors considered only one choice of the strength of IR, in 
the present study we will consider three choices ranging from a 
weak to a fairly strong flux, similar to that emitted by a typical 
young star-cluster. The different strengths of flux are acquired by 
varying the temperature of the source of IR and the rate of photon- 
emission, as recorded in Table 1. We will then compare the evol- 
ution of an irradiated, turbulent cloud with that of a similar cloud 
allowed to evolve only under self-gravity (and no IR). Of particular 
interest is the spatial distribution of dense pockets of gas within the 
MC in each realisation. The plan of the paper is as follows. We shall 
begin by briefly discussing the expected fate of irradiated MCs for 
different choices of gas density. The numerical scheme adopted for 
this work will be described briefly in §3, and the simulations will 
be discussed in §4 before results are formally presented in §5. We 
will conclude in §6. 



2 EVOLUTION OF AN IRRADIATED CLOUD 

The physical situation is illustrated in Fig. 1, where the projection 
(on the plane of the sky) of the cloud facing the external source 
of radiation, S, and located a distance r s from S have been shown. 
The flux of IR, dZ, received by a small area element on the surface 
of the cloud as that shown by the small shaded portion in Fig. 1, 
located at a distance r from the source, S, is 
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where higher order terms involving (R c ulr s ) have been neglected 
as r s » R c id in general. The force, F r , exerted by this flux of IR on 
the cloud is then 
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where c is the speed of light, n,, the average number density of 
ionised gas and £ <2) , the second recombination coefficient defined 

as 
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In the equation above, [} 



Ht^/kT" ' * s me equilibrium tem- 
perature of the ionised gas, and ipiiP) is the second recombination 
coefficient (Spitzer 1978). (Z= 1, for a fully ionised gas). 
Using Eqn. (1), the expression for F r above becomes 

Note that the denominator in this expression for radiation-force has 
units of velocity that may be interpreted as the rate at which the 
ionisation-front advances. As expected, the force becomes vanish- 
ingly small at large distances, r s . Equation (3) for the force exerted 
by the flux of incident radiation permits us to define the momentum, 
p r , delivered by this flux to the gas within the cloud; 

w - dPr 
r ~ dt' 
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For a molecular cloud composed of the usual cosmic mixture the 
average mass of gas particles within it, m =4xl(T 24 g. If the irradi- 
ated cloud loses gas at roughly constant velocity, v, then the rate of 
mass-loss, M/ os ,(r), for this cloud can be obtained by applying the 
conservation of momentum to the mass lost, M/ OM (r), so that 



M,,,„(f) 



dM hss (t) 
dt 
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The expression of the mass-loss rate here assumes that the 
density, n, within the cloud is uniform and that there is no second- 
ary loss of the IR-flux due to absorption by dust within the cloud. 
In a real cloud though, a fraction of the incident flux would be lost 
towards heating the dust within the cloud. The heated dust will of 
course, re-emit in infrared wavebands. The assumption of a pristine 
cloud devoid of any dust, however, is unlikely to alter significantly 
the order-of-magnitude estimate provided by Eqn. (4). The times- 
cale, t, over which an irradiated cloud is likely to survive is simply, 

Mm 

t = - . (5) 

Mfe M (f) 

Using Eqn. (4) we calculate the mass lost by irradiated clouds 
of different average densities, n, over their respective free-fall 
times. The results of a demonstrative calculation performed for a 
photon-ionised gas of average density, n, ~ 10 cm 4 , and main- 
tained at an equilibrium temperature, T ion = 10 4 K, have been plot- 
ted in Fig. 2. The test cloud for each choice of flux-strength was 
kept at a distance r s from the source such that R c id/r s = 0.01. Start- 
ing from the lowest characteristic upward, the average gas density 
within a cloud, n, increases by an order of magnitude. The three 
green vertical lines for the respective choice of Nt y c divides the 
characteristics into survival-regions, i.e., the range of IR-fluxes that 
a cloud can possibly endure. The minimum mass of a cloud that can 
survive without being photoevaporated by the incident radiation- 
flux is defined by the intersection of these vertical lines with the 
characteristics on the plot. 

Small rarefied clouds sit leftward of the vertical line for 
Ni.yc = 10 48 s _I and therefore, are likely to be photo-evaporated 
even by a relatively weak flux of IR. Denser and more massive 
clouds, to the right of this line, will be steadily photo-evaporated 
by a stronger flux. Photoablation is therefore the likely fate of pro- 
gressively larger clouds exposed to a flux of ionising radiation. 
Some of the more massive clouds occupy the top right-hand corner 



of the plot to the right of the vertical line for Ni y c 



10 5 



and evidently need stronger fluxes before they can start losing their 
mass. 

The test cloud used in the simulations below (n ~ 10 4 
cm" 3 ,/? c « = 1 pc; see §3.4 below), lies on the characteristic just 
below the uppermost in the region between the fluxes correspond- 
ing to 10 49 s _I and 10 51 s _I , and therefore, is expected to steadily 
photo-ablate. In fact, the mass-loss rate for this test cloud according 
to Eqn.(4) above is a few times 10~ 4 M G yr" 1 for N LyC ~ 5 x 10 49 
s , and R c id/r s ~ 0.08. The mass-loss rate for this cloud suggests, 
it is likely to survive at least an order of magnitude longer than its 
free-fall time. 



3 NUMERICAL METHOD 

3.1 Smoothed Particle Hydrodynamics (SPH) 

The SPH algorithm was first introduced by Gingold & Monaghan 
(1977), and Lucy (1977) to handle complex problems in astro- 
physical context and other areas of computational fluid dynamics 



MC 




Figure 1. A projection showing the location of the molecular cloud (MC), 
of radius R c u, relative to the source of ionisation, 5 . The amount of flux 
received by the shaded area element on the surface of the cloud is given by 
Eqn. (1) in the main text. 




9 9.5 
log 10 (Flux/[cm" 2 s" 1 ]) 

Figure 2. A simple classification scheme for MCs that enables us to predict 
its evolution when exposed to a flux of IR; the calculation assumes «, = 10 
cm -3 , Tj on = 10 4 K and r, / R c /d=l2. The three choices of NtyC used in the 
present work are shown by green vertical lines. The intersection of these 
lines with the cloud mass-flux characteristics defines the minimum mass 
required for a cloud of given density to survive the incident flux of IR. See 
text for description. 



(e.g. Monaghan 2005). We use our numerical hydrodynamics code, 
SEREN, an open-MP parallised algorithm that has been extensively 
tested for numerous applications of astrophysical interest (Hubber 
et al. 201 1). The fundamental quantity in an SPH code, the density 
of particles, is calculated using the gather technique. The search 
for nearest neighbours of an SPH particle, and the calculation of 
the net force (gravity+hydrodynamic) on an SPH particle is done 
by distributing particles on the Barnes-Hut tree, with a cell-opening 
angle of ~0.2; quadrupole moments of distant cells in the tree struc- 
ture are also included in the calculation of forces. We have used the 
standard prescription for SPH artificial viscosity with the viscous 
parameters having values (0.1,0.2). 

The effects of ionising radiation are included through a ray- 
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Table 1. Listed below are the choices of ionising flux and the temper- 
ature of the source used in various test cases (Spitzer 1978). The tem- 
perature of the ionised gas is maintained at, T lo „ = 8000 K in cases 2, 3 
and 4. 



Serial 
No. 




r -2 -Ii 

[cm s J 


1 


No ionising radiation 




2 


4xl0 48 ; 
36,000 K 


1.10 X 10 7 


3 


5xl0 49 ; 
47,000 K 


1.37 X 10 s 


4 


7xl0 51 ; 
51,000 K 


1.93 X 10 10 



A/iyCi T s1ar , "£-Nlvc < aTe respectively the rate of ionising photons, the 
temperature of the source star and the flux of ionising radiation 

casting algorithm devised by Bisbas et al. (2009), this scheme is an 
inherent part of SEREN. The basic features of this algorithm are : it 
- (a) determines the position of the ionisation front (IF) at a given 
instant of time, and (b) assigns appropriate temperatures to the SPH 
particles. Of the two, (a) is achieved via the HEALPix algorithm 
(Gorski et al. 2005). Individual rays originating from the source 
of IR that trace the direction of propagation of the radiation are 
distributed uniformly on a sphere by HEALPix. Rays are distributed 
at different levels, /, between and 7 so that there are 12 X 4' rays 
on any level. The angular resolution of the IF is determined by the 
value of /. The maximum intensity, I max , at the IF from Eqn. (6) is 

_ m 2 N Ly c 

Imax = 4^(2) ' (6) 

where m = 5, is the mean mass of a gas molecule and X=0.7, 
and £ {T> is the second recombination coefficient defined by Eqn. (2) 
above. 

The integral form of the above expression for intensity at a 
point r is, 

l(r) = Jp 2 (re)r 2 dr, (7) 

where e is the unit vector in the direction of r. Equation (7) is solved 
iteratively to determine r such that I(r) = /„ rav , while calculating the 
SPH density at each point of the iteration. 



3.2 Thermodynamics 

All points satisfying the inequality I(r) < I max , lie within the HII 
region and are therefore assigned a temperature, T IO „, the equilib- 
rium temperature of the HII region. While the equilibrium temper- 
ature probably depends on the coupling between matter and radi- 
ation (e.g., Genel et al. 2012), in the present study, T lo „, is assigned 
a fiducial value of 8000 K. Neutral gas within the irradiated cloud 
is treated isothermally and maintained at a temperature, T„ = 20 K. 
The algorithm also includes a temperature smoothing scheme that 
alleviates the discontinuity arising due to a sharp difference in the 
temperature of the ionised and neutral gas particles. The smooth- 
ing scheme is essentially a Taylor-series approximation to the first 
order. 



3.3 Sink particles 

An SPH particle with density higher than the threshold, pthresh ~ 
10" 15 g cirr 3 , is replaced by a special type of particle, the sink, 
following the prescriptions made by Bate & Burkert (1997). Apart 
from the density threshold, a sink particle is also has a radius which 
in the present simulations is set as 2.5 times the SPH smoothing 
length of a prospective sink particle. SPH particles bound gravita- 
tionally to the sink are accreted by it. The sink particle in the present 
numerical exercises represent a protostellar object. The minimum 

mass of an SPH particle, M„„„ = N^a/ where N,,^ = 50 
and N gaj = 550, 000, are respectively the number of nearest neigh- 
bours of an SPH particle and the number of particles within the 
cloud. With this choice of SPH parameters the ratio of the thermal 
Jeans mass, Mj(p ~ 10~' 6 g cirr 3 , T„ = 20 K) against the min- 
imum mass, M,,,;,,, is 6 which satisfies the Bate-Burkert criterion of 
resolving the gravitational instability. 

3.4 Initial conditions 

The test cloud used in the present study was modelled as a uniform- 
density sphere of unit mass and radius and was assembled by ran- 
domly positioning particles within it. Since the randomness in pos- 
itioning particles is inherently associated with Poissonian noise, we 
first evolved it for a fraction of a sound-crossing time allowing the 
numerical viscosity to dissipate the spurious noise. This procedure 
settled the initial assembly of particles which was then rescaled to 
the desired dimensions. The mass and radius of the cloud are re- 
spectively, M c id = 400 M G , and radius, R f u = 1 pc, where as the 
gas within it was maintained at a uniform temperature, T cH = 20 K. 
A supersonic, Gaussian velocity field with random amplitudes and 
a relatively steep power-spectrum, P(k) oc kr 4 ', was then superposed 
on this cloud Q This velocity field had an initial Mach number of 
10 and was first set up on a grid with 128 3 cells in the fc-space. The 
test cloud was then mapped on the i-space and velocity compon- 
ents for individual particles were calculated via interpolation. Then 
a random phase was added to each velocity component of a particle 
while transforming it back into the Cartesian-space. The cloud with 
its turbulent velocity field was then allowed to evolve under self- 
gravity; this is the first case listed in Table 1. This cloud soon de- 
veloped filamentary structure as shown in the rendered density plot 
of Fig. 3 {t ~ 0.05 Myr), which was then used as the initial condi- 
tion in the remaining 3 cases, listed 2 through to 4 in Table 1 . The 
cloud age was not reset to zero at the time of introducing the source 
of IR in these latter cases. Each test case discussed here took a little 
over 5000 CPU hours and was run on the HYDRA supercomputing 
cluster at the Indian Institute of Astrophysics. The cluster is com- 
posed of the Intel Xeon-5675 processor. Simulations were allowed 
to form 150 sink particles before they were terminated. 



4 NUMERICAL SIMULATIONS 

4.1 Turbulent cloud without a source of ionising radiation 

Collisions between random gas flows within a turbulent cloud dis- 
sipate turbulent energy and generate filamentary structure within 
the cloud as can be seen in the renders plotted in Figs. 3 and 4. The 
initial velocity field, with a power- spectrum considerably steeper 

1 k&P(k), are respectively the wave-vector and the power in the wave- 
vector space 
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Figure 3. A rendered density plot, at t = O.ltff ~ 0.05 Myr, showing the 
mid-plane of an initially turbulent MC. Injected turbulence readily gener- 
ates dense filaments within the cloud, as can be seen in this plot. This cloud 
is then used as the initial condition for the remaining test cases when the 
source of IR is turned on (ter ~ 0.27 Myr). 

than the Kolmogorov type, was injected between wavenumbers 1 
and 8. As a result turbulent energy cascades down via shocks from 
larger spatial scales. Respective plots in the panels of Fig. 4 show 
the rendered density images of the mid-plane of the cloud at suc- 
cessive epochs. After a little over one free-fall time the original 
cloud collapses to a centrally located dense filament with a few ar- 
terial extensions. The advanced stage of the collapsed cloud can be 
seen more closely in the rendered plot of Fig. 5 where positions of 
a few protostellar objects have also been marked. As expected, the 
formation of protostars proceeds along the dense filaments. 

4.2 Irradiated turbulent cloud 

Exposure to a flux of ionising radiation commonly causes the mo- 
lecular cloud to lose mass via photo-ablation. Unlike small, dense 
clumps which when exposed to IR implode to form comet-shaped 
globules, the larger clouds such as the one considered in the present 
study develops sporadic pockets of dense gas, often filamentary and 
eventually spawn stars. In this next set of 3 simulations we endeav- 
our to study the effect of varying the strength of IR on the evolution 
of a MC. 

4.2.1 Weak ionising radiation 

The exposure of a MC to a relatively weak flux of IR, listed as case 
2 in Table 1, produces a C-shaped ionisation front (IF), on the sur- 
face of the cloud exposed to the source of IR. The incident flux of 
IR also shocks the pre-existing filaments within the cloud, visible in 
the plots on the upper two panels of Fig. 6. Interestingly, the shock- 
ing of a clump in one of these filaments also generates a small, 
faint elephant trunk-like extension that can be seen in the lower 
left-hand corner of the second panel of this figure. The growth of 
this trunk is probably the result of a backward re-expansion of ma- 
terial from a shocked clump in a filament within the cloud. This 
possibility has previously been demonstrated in numerical simula- 
tions by for instance, Williams et a/.(1999), and more recently by 




Figure 4. A sequel of column density plots showing the mid-plane of the 
turbulent cloud at different epochs; t = 0.1,0.15, 0.2 & 0.25 Myr, for images 
in respective panels from the top to bottom (tf/ ~ 0.27 Myr). Turbulence 
within the cloud generates filamentary structure that can be seen in the fig- 
ure on each panel. 
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Figure 5. A column density plot showing an advanced stage in the evolution 
of the turbulent cloud in case 1; t = 0.31 Myr~ 1.2 f//. The position of the 
sink particles within the collapsed cloud has been marked with * on this 
plot. Note that sinks represent protostellar objects. 

Gritschneder et al. (2009,2010). Protruding trunks are commonly 
found in star-forming clouds and some even harbour a YSO in the 
head region such as those reported in the famous M16 (Eagle) neb- 
ula (e.g., White et al. 1999). However, the trunk observed here ab- 
lated soon after its appearance. Turbulence injected by the flux of 
IR further assists the formation of filaments within the cloud, the 
late phase of which, with a network of filaments and star-forming 
sites within it can be seen in the rendered density plot in Fig. 7. 
Also marked on this rendered plot are the positions of protostars 
that appear aligned with dense gas filaments. 

The fraction, / = M re ,(t)/M c ij, of the irradiated cloud retained 
has been plotted in the upper panel of Fig. 8. The fraction of the 
cloud lost is therefore (1 - /(?))■ Exposure to a weak flux of IR 
causes the cloud in this case to lose ~30% of its original mass and 
so a large portion of the cloud is still retained. The residual cloud 
evolves like the turbulent cloud discussed in §4. 1 above, although 
the network of filaments within it, as can be seen in Fig. 7 and 
the lowest panel of Fig. 6, is significantly more branched in stark 
contrast to the filament seen in Fig. 5. The incident flux of a weak IR 
appears to drive many more modes. The simulation was terminated 
after about 0.9 free-fall times (t = 0.24 Myr), by which time 150 
protostars had formed within the cloud. 



Figure 6. The mid-plane of the cloud in case 2 irradiated by a weak flux 
of IR is shown in this sequel of column density plots. A faint trunk that 
eventually evaporates due to the flux of IR can be seen in the lower left 
corner of the second panel. While exposure of the cloud to the incident flux 
of IR causes the cloud to loose mass, the radiation-induced shocks produce 
dense filamentary structure within the cloud; especially see the plot in the 
lower panel (t = 0.216 Myr -0.8 tff). 



4.2.2 A relatively strong ionising radiation 

With a stronger flux of IR, the MC in this case, listed 3 in Table 1, 
loses mass at a much higher rate than that observed in the previous 
realisation. Consequently, after t = 0.8?// ~ 216 kyrs, a little over 
40 % of the original cloud is lost via photo-ablation. Filaments in 
the original cloud, as is evident from the sequel of rendered dens- 
ity plots in Fig. 9, are shocked by the incident flux of IR. This flux 
also produces a thin, dense shocked-shell on the surface of the ir- 
radiated cloud. The surface of the IF as can be easily seen, is cor- 
rugated. A closer examination reveals its wiggly nature which re- 
sembles features of the well-known thin-shell instability (TSI). The 



TSI, as the name suggests, is triggered in thin shells confined by 
ram-pressure either on one or both faces. The instability grows due 
to the transfer of momentum between perturbed regions of the shell 
as demonstrated by Vishniac (1983), through a detailed perturba- 
tion analysis of the problem. The dynamically unstable nature of 
the irradiated surface has also been discussed by Williams (2002), 
who demonstrated the appearance of corrugations on the surface of 
the IF. However, it is unclear if these later corrugations are also the 
manifestations of the TSI, though for the moment we do not attempt 
to distinguish between the two. 

We therefore use a previously derived expression for the 
length of the fastest growing mode, Af ast , to verify if the instabil- 
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Figure 7. A column density plot showing the picture in the lowest panel 
of Fig. 6. Dense filaments within the irradiated cloud are readily visible 
in this plot and a few sink particles (i.e., protostellar objects) have been 
marked with *, which unsurprisingly, are located along the filaments (t ~ 
0.24 Myr). 



.=3.6e4 K 
-4.7e4 K 
,=5.1 84 K ■ 



ity is indeed resolved in this simulation. The length of the fastest 
growing mode in a shocked layer of uniform surface density, cr , 
calculated by Anathpindika (2010) is 



l fast 



In 

kfast 



where 



K fast 



ttGctq 



This mode grows on a timescale, t 



(8) 



growth 



A fast I a; a 



(ksT^ Jihu) 112 is the sound speed and, T gas , the average temper- 
ature of the gas in this layer, which is ~ 10 4 K. The radius of the 
irradiated cloud, R' 2 = R M + dR, where R cld is the radius of the 
original cloud and dR is the thickness of the shell. In the present 
case, dR/R' cld ~ 10~ 2 , and PeIPs ~ a\ u is the square of the sound- 
speed in ionised gas. Plugging in the appropriate values we get, 
A f as , ~0.05 pc, and t groMn ~0.22 Myr, which agrees with the times- 
cale on which wiggles first appear on the surface of the shocked- 
shell. 

The ratio of Af ast to the average SPH smoothing length, X, is 
a good indicator of the spatial resolution, which in the present case 
is 3. It is therefore likely that the instability has been only barely 
resolved here. Better resolution of the instability demands at least 
an order of magnitude increase in X, which transforms in to about 3 
orders of magnitude increase in the number of SPH particles. This is 
computationally demanding. However, since we are not concerned 
about the detailed modelling of this instability, the present choice 
of spatial resolution is sufficient. It is therefore clear that radiation- 
induced shock produces pockets of dense gas along the irradiated 
surface, and further, a few other isolated pockets of dense gas can 
also be seen along the shocked filaments in the interior of the cloud. 
The rendered density plot in the left-hand panel of Fig. 10 shows 
the irradiated cloud more closely along with the location of sink 
particles on it. As with the earlier case, this simulation was also ter- 
minated at / ~ 0.22 Myr, by which time the cloud had formed 150 
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Figure 8. Top-panel : Plot showing the fractional mass, /, of the irradiated 
cloud retained in each of the three realisations. Lower-panel : The mass of 
gas within the cloud above the density threshold, Pthresh ~ 10~ 18 g cm~ 3 , 
has been shown at different epochs of time. Relatively stronger shocks, such 
as the ones in cases 3 and 4, assemble material in the dense phase on a much 
shorter timescale. 



sink particles. We note that the sink-formation timescale in cases 2 
and 3 is therefore mutually comparable. 



4.2.3 Strong ionising radiation 

In this final test case, listed number 4 in Table 1, the flux is the 
strongest with source characteristics similar to those of a typical 
young star-cluster. The irradiated cloud in this case evolves in a 
manner similar to that observed in the previous case, albeit on a 
much shorter timescale and has a much higher mass-loss rate, as 
also reflected by the plots in Fig. 8. Analogous to the previous 2 
cases, exposure to a flux of IR shocks the surface and the interiors 
of the cloud. Some of the clumpy features near the surface briefly 
appear to develop trunk-like features that are stymied soon due to 
rapid photo-ablation of the surface. The surface of the irradiated 
cloud, as in case 3, also shows evidence of the thin-shell instabil- 
ity discussed previously in §4.2.2. Since the broad features of the 
irradiated cloud in this case match with those in the previous real- 
isation, we only show the final state of the cloud in the right-hand 
panel of Fig. 10. Unlike the irradiated cloud in the previous case 
the initial star-formation in this case is confined predominantly to 
the shocked shell 
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Figure 9. Mid-plane of the irradiated cloud in case 3 has been shown in these rendered column density plots. Exposure to a relatively strong flux of IR ablates 



the cloud while also shocking the dense filaments within it as can be seen in the plots on each panel; 
the left-, central and the right-hand panels respectively. 



t ~ 10 3 , 2xW and 2.1xl0 3 yrs (~ 0.7?//), for plots in 




Figure 10. Left-panel: A rendered column-density plot showing the mid-plane of the irradiated cloud in case 3 (/ ~ 0.75f//), after the formation of the first 
few protostellar objects. The protostars, represented by the sink-particles, have been marked with '*'. Observe that the formation of stars commences in the 
shocked layer on the cloud surface, though another potential site of star-formation appears in the central region of the irradiated cloud. Right-panel: A rendered 
column density plot showing the mid-plane of the irradiated cloud in case 4. As in case 3, most of the star-formation in this case is also confined to the shocked 
layer on the irradiated surface of the cloud (t ~ OAMyr = 0.33?//). Shocked filaments within the cloud, thus far, do not appear to have acquired sufficient 
density to spawn stars. 



5 DISCUSSION 

Star-formation according to Krumholz & McKee (2008), is un- 
likely to commence unless gas density is higher than a threshold. 
The proposition is attractive, for such a threshold will essentially 
reduce the problem of star-formation to simply an investigation 
into the processes likely to assist the assembly of a sufficiently 
dense volume of gas. However, it is common knowledge that the 
dynamical state of gas is also inextricably interwoven with its tem- 
perature, where for convenience, we include turbulence as just an- 
other thermal component. It is therefore difficult to define a density 
threshold for the formation of stars. 

There are two crucial issues before us : (i) the process, usually 
referred to as a trigger, that assembles pockets of dense gas, and (ii) 
the processes that maintain these pockets at the right temperature 
so they can spawn stars. Clearly, this latter question is beyond the 
remit of the present work, however, by maintaining isothermality 
within the cloud, at least the first half of the problem can be ex- 



amined. Other crucial points that are also discussed while explain- 
ing the formation of stars are of course, the timescale over which 
stars form which is also an indicator of their age, and the distribu- 
tion of stellar masses at their birth, characterised by the initial mass 
function (IMF). 

We now examine the results from our test simulations to find 
answers to these questions. While the turbulent velocity field ori- 
ginally introduced into the cloud produces dense filaments of gas, 
the effect of the incident flux of IR on the cloud depends on its 
strength as we have discussed so far. Firstly, all the 4 test simula- 
tions show that a suitable triggering mechanism enhances the rate 
at which gas is assembled in dense pockets. While some pockets of 
dense gas are isolated, yet others lie in larger filamentary structures. 
The strength of the trigger, however, determines the essential phys- 
ical characteristics such as the mass of these pockets of dense gas 
and the timescale over which they are assembled within the cloud. 

The plots in Fig. 8 thus offer an instructive comparison. The 
retention factor, /(M), plotted for the irradiated cloud in each test 
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Figure 11. Upper-panel : Figure showing the probability distribution func- 
tion (PDF) for gas density within the cloud at an advanced stage of its evol- 
ution in each of the 4 test cases. A log-normal distribution fits all the PDFs 
reasonably well, though they also have a prominent power-law tail. We note 
that the width of the log-normal distribution is narrower for the irradiated 
clouds and the narrowest for the choice of the strongest flux. The best- fitting 
log-normal distributions have also been shown for comparative purposes. 
Lower-panel : The power-spectrum for each of the 4 test cases has been 
plotted. The spectra show remarkable similarity to the Kolmogorov spec- 
trum. We note that maximum power is cascaded down to smaller spatial 
scales in case 4. 



case shows that a progressively smaller mass of the cloud is re- 
tained with increasing strengths of the incident flux. The observed 
rate of photoablation in cases 3 and 4 is a few times ~ 10~ 4 M Q 
yr~', which is consistent with that calculated using Eqn. (4) above. 
Plotted in the lower-panel is the timescale on which gas within the 
cloud overhauls a fiducial density threshold, p t hresh ~ 10 -18 g cm -3 . 
Thus, with relatively weaker shocks resulting from damped turbu- 
lence in case 1, it takes considerably longer, more than 70 % of the 
free-fall time, to overhaul a fiducial density threshold of 10" 18 g 
cm 4 ; see the plots in the lower-panel of Fig. 8. 

Exposure to a weak flux of IR as in test case 2 makes little dif- 
ference to the rate at which gas overwhelms this threshold. Increas- 
ing the strength of the IR though strongly ablates the cloud and, 
more importantly, gas within the cloud also enters the dense phase 
on a relatively shorter timescale. This can be seen easily from the 
respective plots for cases 3 and 4 in the lower-panel of Fig. 8. The 



mass of gas within the irradiated cloud above the density threshold, 
M(p,;, rf ,,;,), rises most rapidly in case 4, where the IR-flux is the 
strongest. 

Our conclusions from the set of simulations discussed above 
are also corroborated by those reported recently by Walch et al. 
(2012). These authors have drawn attention to the dependence on 
the fractal dimension of the evolution of an irradiated MC. They 
suggest two possible scenarios for a fractal cloud. Thus, a MC with 
a relatively large fractal dimension (>2.6) was observed to pro- 
duce pillar like features where as the cloud with a fractal dimension 
smaller than ~2.2 became perforated as dense structure within the 
cloud was shocked by the incident radiation. To this end our res- 
ults are consistent with this latter scenario; see the rendered plots 
in Figs. 6 and 9. 

In another related work Dale & Bonnell (2012), while report- 
ing the findings from their simulations contradict these results, in- 
cluding ours. They found that the evolution of a turbulent MC, a 
few thousand M Q s massive, was largely unaffected even after ex- 
posing it to a strong flux of IR. Thus, irrespective of the flux of IR, 
the turbulent MC in their study ended up as a single central fila- 
ment, although trunk-like features developed on the surface of this 
filament when exposed to IR. However, the resulting MFs were mu- 
tually comparable. Dale & Bonnell (2012) therefore suggest that 
even a relatively strong flux of IR is unlikely to significantly alter 
the evolution and the star-formation history of that cloud. 

This radically different result is probably due to - (a) the dif- 
ference in the modes in which turbulent energy was initially injec- 
ted, and (b) a relatively larger initial mass and density of their test 
cloud that could have impeded radiation-induced turbulence. The 
formation of a single filament in their simulations also suggests 
that the mode driven by the initial turbulence was dominant and 
grew rapidly before the energy injected by the IR-flux could have 
any impact. We now discuss the principle findings of our study. 



5.1 Probability density function (PDF) 

The probability density function is a useful diagnostic to study 
the effects of physical conditions on the distribution of gas dens- 
ity within a MC. The PDFs for each of the 4 test cases have been 
plotted in the upper-panel of Fig. 1 1 . In general, the PDFs appear 
to be a combination of power-law and lognormal distributions. It is 
evident that gas distribution within the irradiated cloud in cases 2 to 
4 spans a wider rage of density as compared to the turbulent cloud 
in case 1 where it was not exposed to IR. The PDF in this latter case 
has distinct power-law tails at both, the low as well as high-density 
ends, albeit the power-law for higher densities is much steeper. The 
distribution for intermediate densities though, is roughly lognormal 
as can be seen in the plot in the upper-panel of Fig. 1 1 . The PDFs for 
the irradiated clouds in the remaining 3 cases are mutually similar, 
however, the power-law tail for each one of them is significantly 
steeper than that observed in case 1. The PDFs derived here and 
those reported by us for a shocked-cloud in an earlier work (Anath- 
pindika & Bhatt 2011), bear striking resemblance with those de- 
duced by Schneider et. al. (2012), for the Rosette MC irradiated by 
an adjacent cluster of stars, NGC 2244. These latter authors along 
with a few others (e.g., Hill et. al. 2011), also suggest the possib- 
ility of the occurrence of more complex PDFs than those derived 
here. 
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5.2 Power-spectrum 

The spectra plotted for each of the four test cases and shown in the 
lower-panel of Fig. 11 resemble the Kolmogorov spectrum. The 
derived power-spectra are also consistent with those for MCs in the 
Galactic ring (Roman-Duval et al. 2011). Progressively increasing 
the strength of the turbulence both, drives a broader range of wave- 
numbers and produces a considerably steeper spectrum such as the 
one for case 4. Turbulent energy in these latter cases of radiation- 
driven turbulence cascades down to smaller length-scales which 
probably induces fragmentation on smaller spatial scales. The ef- 
fect of the energy-cascade also manifests itself in the PDFs plot- 
ted in the upper-panel of Fig. 11 and discussed in the previous 
section. Enhanced fragmentation also shifts the protostellar mass- 
distribution towards lower masses as is indeed observed in the mass 
distributions derived here and discussed in §5.4. 



5.3 The star-formation rate: Is triggered star-formation 
rapid ? 

Having discussed one aspect of the problem, we now move to the 
next where the possible effects of shocks on the distribution of stel- 
lar masses are examined. The accretion history of sink particles 
in each simulation is plotted in Figure 12, where the combined 
mass of sink particles in each realisation has been shown. The 
spread along the abscissa reflects the distribution of ages of the 
sink particles in a particular test case. Importantly, the steep nature 
of the plots suggest that the formation of the first sink particle is fol- 
lowed by an episode of rapid sink-formation and a number of them 
form over the next few 10 4 years. Each simulation was terminated 
after the formation of 150 sinks though this means, the simulations 
were not coeval at the time of termination. We observe that sink- 
formation timescale in each test case is only a fraction of the age of 
the cloud itself, or indeed, its freefall time. 

The rate at which ordinary gas is converted into stars is quan- 
tified by the so called, star-formation rate (SFR). The SFR is also a 
measure of the efficiency with which an external trigger assembles 
potential star-forming sites. The star-formation rate (SFR), we re- 
mind, is actually the rate at which sinks form and has been calcu- 
lated for each test case; see Table 2. It is the largest for case 1 where 
the turbulent cloud was left by itself to evolve under self-gravity. 
Irradiated clouds in cases 2, 3 and 4, on the other hand, show a 
mutually comparable, though relatively low SFR. However, these 
latter values of the SFR are consistent with those derived for typical 
star-forming clouds in the nearby universe (e.g. Lada et al. 2012). 
The relatively small SFR in cases 2, 3 and 4 is not entirely sur- 
prising since injected turbulence causes sub-fragmentation of lar- 
ger structures to produce smaller cores and therefore, less massive 
stars. This tendency is further reflected by the distribution of sink 
masses plotted in Fig. 13. 

Authors have often cited similar values of the SFR to support 
the hypothesis of slow star-formation (e.g., Krumholz & Tan 2007). 
Conversely, we have shown here that shocks triggered due to con- 
tinuously driven turbulence rapidly assemble gas in the dense phase 
before eventually spawning stars. Stars thus form on a relatively 
short timescale, only of the order of a fraction of the free-fall time 
of the original cloud, though protostellar feedback may further re- 
duce the SFR and so, the derived values here are likely tobe the 
upper limits. 

Recent studies of a number of star-formation sites such as 
those within the Rosette nebula (e.g., Wang et al. 2008, Schneider 
et al. 2012), the RCW34 (Bik et al. 2010), or the SF038 (Choud- 



Table 2. The observed rate of sink-formation for each test case (SFR = 
Z M„" h ); t $ i n k s marks the epoch when a simulation was terminated, in 
other words the age of the cloud in each case. 
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SFR[xl(T 4 ] 




[M Q yr-'] 
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5.61 
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1.07 
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1.42 


4 


1.94 
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Figure 12. Plots showing the rate at which the total sink-mass in each of 
the 4 test cases increases. Observe that sink-formation in irradiated clouds 
occurs on a relatively short timescale; it is, however, significantly delayed in 
the cloud not exposed to the IR-flux where sinks appear only after a free-fall 
time (t ff ~ 0.27 Myr). 

hary et al. 2010), have revealed several interesting properties, some 
of which can be reconciled with the numerical models discussed in 
this study. Yet in the well-known MCs such as the Orion and Mono- 
ceros, Wilson et al. (2005) have reported a systematic increase in 
the age of stellar populations, indicative of sequential propagation 
of star-formation. Unfortunately, due to the absence of a feedback 
mechanism in our models we are ill-qualified to comment on this 
aspect of the star-formation process. At the least, we add that se- 
quential propagation of star-formation is likely via feedback from 
the first generation of stars within a star-forming cloud. 

5.4 Mass function (MF) 

The mass function(MF) generated using 150 protostellar objects 
(i.e., sinks) for each test case is shown in Fig. 13. It will be useful 
to recall that the age of the fragmented cloud, and therefore the MF, 
at this stage is different in each test case, varying between ~0.33 
tff in case 4, to —1.2 tff in case 1. This can also be seen from the 
sink-accretion histories plotted in Fig. 12. Quite clearly, the MF in 
each test case has a power-law nature of the type 



for a non-zero a. The parameters of the MFs derived for each case 
tested in this study are listed in Table 3. The MF for each case 
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Figure 13. Top panel: The sink mass function(MF) for cases 1 and 2 are shown in the left and right-hand plots. Lower panel: Plots in the left and right-hand 
panels show the sink-mass functions for cases 3 and 4. Each MF was generated using masses for 150 sink particles. Note that MFs in the first 3 cases are 
incomplete below ~0.05 M Q , where as in case 4, the MF turns-over towards masses smaller than ~0.2 M Q . 



was derived by convolving 10 s realisations of the respective data 
set within \a. Starting with a relatively shallow power-law in case 
1 where turbulence is continuously damped, the MF steepens pro- 
gressively as turbulent energy is continuously injected by increas- 
ingly stronger radiation-driven shocks in the remaining 3 cases. 

The MF in case 4 where the flux of IR is the strongest, at- 
tains the canonical slope at the high-mass end. It develops a knee 
at ~0.4 M Q , before turning over into lower masses at ~0.2 M Q . The 
slopes, (a i, 0-2), of this MF for masses <0.4 lie within the error-bars 
on the stellar initial mass function (IMF) due to Kroupa (2002). 
On the other hand, in cases 2 and 3 where the flux of IR is rel- 
atively weaker, the MFs develop a knee at ~0.2 M Q . However, it 
shifts right-ward, i.e., to a slightly higher mass, ~ 0.4 M Q , for case 
1 where the derived MF itself is shallow and the cloud is not ex- 
posed to radiation. We note that the minimum mass, M„„„, of an 
SPH particle defined in §3.3 above sets the threshold for the com- 
pleteness of the MFs derived here. In the simulations discussed in 
the present work, M„„„ ~0.04 M , rendering the MFs incomplete 
below ~0.05 M . The minimum mass, M„„„, we remind, has been 

defined in §3.3 above as N Be «,s( It is therefore not possible for 
us to comment on the possibility of the formation of stars smaller 
than ~0.05 M via turbulent fragmentation. 



6 CONCLUSIONS 

In this article we presented a comparative study of the evolution of 
turbulent MCs. Turbulence, irrespective of whether it is damped or 



injected continuously, generates filamentary structure within MCs. 
In the first test case of this study, turbulence was allowed to damp 
without ever being replenished where as in the remaining three 
test cases, it was injected continuously via radiation-driven shocks. 
These latter three cases show that irradiating a large MC profoundly 
affects its evolution, and therefore, the star-formation history. The 
cloud with damped turbulence, however, evolves significantly dif- 
ferently compared to that in which turbulence is continuously injec- 
ted via an external source. Not only is there a difference in the evol- 
ution, but the timescale on which the cloud evolves also varies and 
in fact, the cloud with stronger turbulence within it evolves much 
faster. Thus in the present work, the injected turbulence within the 
cloud is the strongest in case 4 and so, star-formation within the 
cloud in this case commences on the shortest timescale in compar- 
ison to the other 3 cases. 

Global gravity soon dominates gas dynamics within the cloud 
evolving purely under self-gravity so that few modes other than the 
Jeans unstable mode are driven. Consequently the Jeans mode be- 
comes the principal mode of instability so that gas within the MC 
is primarily channelled into this mode. Furthermore, this tendency 
is reflected by the PDF as well as the MF for the MC in case 1. 
Without much surprise, the PDF for this MC is largely predisposed 
towards somewhat higher masses, an obvious mis-fit to the canon- 
ical IMF. Driven turbulence, on the other hand, injects energy in to 
various unstable modes which aids formation of dense structure. In 
this study turbulence was injected within MCs by an uninterrupted 
source of IR. Simulations were performed by varying the strength 
of the IR-flux, spanning a range of radiation emitted by a typical 
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Table 3. Physical parameters for the mass function (MF) derived in each 
test case. The indices a[ and 0-2 are the respective slopes on either side 
of the knee (i.e., Mj/nt < Mk nee and M.„-„^ > Mknee) of the MF. The MF 
only in case 4 turnsover for (M„„t/Mo) < 0.2, and the corresponding 
slope is denoted by 03 . 



Case 


<*l 


(i'2 


knee 
[Mo] 


<*3 


turn-over 
[M Q ] 


1 


0.9 


1.6 


0.4 


MF incomplete 
below -0.2 M 




2 


1.3 


2.0 


0.2 


as above 




3 


1.6 


2.2 


0.2 


as above 




4 


2.0 


2.4 


0.4 


-1.3 


0.2 



B-type star to that emitted by young OB-associations. There is little 
doubt that progressively stronger radiation drives shocks of ever 
increasing strength. 

A common feature of the irradiated MCs observed in each of 
the three test cases, 2, 3 and 4 of this study, is the appearance of a 
shocked shell on the irradiated surface facing the source of IR and 
the shocked internal filamentary structure generated by the initial 
turbulence. The secondary shocks driven by the flux of IR create 
new sites of star-formation, particularly in the shocked shell and 
junctions of filaments within the MC. Evidently shocks distribute 
gas over a wider range of density as reflected by the correspond- 
ing PDFs. The range over which gas-density is distributed becomes 
wider with increasing strength of the radiation-driven shocks. In- 
cidentally, these PDFs are consistent with those reported recently 
for star-forming clouds exposed to strong fluxes of IR. 

Multiple stellar-systems appearing in isolated pockets of gas 
compete for gas as envisaged in the standard competitive-accretion 
scenario. The protostellar MFs derived in the present work, though 
fitted by a power-law, are not coeval. Interestingly, the MFs for the 
irradiated cloud in cases 2, 3 and 4 better agree with the canonical 
IMF, or even its more recent variation suggested by Kroupa. In fact, 
in case 4 where the radiation-induced shocks were the strongest, the 
MF even demonstrated a turnover at ~0.2 M which is consistent 
with the Jeans mass for the shocked gas maintained at its pre-shock 
temperature. The MF in case 1, though a power-law, is consider- 
ably shallow.Thus we have shown that cloud models with injected 
turbulence can possibly reconcile some key dynamical features of 
typical star-forming regions. This study, for instance, also shows 
that star-formation once triggered, proceeds rapidly and newer stars 
form within just over a few 10 4 years, which is significantly smal- 
ler than the free-fall time of a typical star-forming cloud. Also, the 
star-formation rate (SFR) derived for the irradiated cloud in cases 
2, 3 and 4 is consistent with that reported for typical star-forming 
clouds. Thus, despite the dissipative nature of turbulence, it can 
profoundly affect the star-formation history of a MC and external 
sources such as a flux of IR can continuously replenish the turbulent 
energy within that cloud. Consequently, principle sources of turbu- 
lence such as proto-/stellar feedback in nearby star-forming clouds, 
or nuclear star-formation activity in galaxies must be poignant for 
controlling the global star-formation history. 
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